General elastic interaction in nematic liquid crystals colloids 
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The new free energy functional that describes general elastic interaction between colloidal particles 
and nematic liquid crystal has been proposed. It generalizes results of the paper [?J on the case 
of arbitrary orientation of colloidal particles and is valid for arbitrary surface anchoring strength. 
Formal analogies and differences between electric particles and colloidal particles in LC are found. 
0^ \ It is first time shown that spur of the quadrupole moment tensor is different from zero SpQ^ 7^ 

in the case of axial symmetry breaking of the director field around colloidal particle in contrast 
to electrostatics. New general presentation for the elastic interaction between arbitrary colloids is 
obtained which contains new multipole terms connected with SpQ^ 7^ that are absent in the 
standard electrostatics. 
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Colloidal particles in liquid crystals (LC) have attracted a great research interest during the last years. Anisotropic 
properties of the host fluid - liquid crystal give rise to a new class of colloidal anisotropic interactions that never 
occurs in isotropic hosts. Liquid crystal colloidal system have much recent interest as models for diverse phenomena 
in condense matter physics. The anisotropic interactions result in different structures of colloidal particles such as 
linear chains in inverted nematic emulsions [l|, [2| , 2D crystals [3j and 2D hexagonal structures at nematic-air interface 

Study of anisotropic colloidal interactions has been made both experimentally 0-0 and theoretically The 
first theoretical approach was developed in [l|, 0] with help of ansatz functions for the director and using multiple 
expansion in the far field area. Authors investigated spherical particles with hyperbolic hedgehog and found dipole and 
quadrupole elastic interactions between such particles. Another approach jg, {| gave possibility to find approximate 
solutions in terms of geometrical shape of particles for the case of small anchoring strength and has provided the 
to connect the type of the interaction potential with the local symmetry of the director field around particles 
The concept of coat has been introduced that contains all topological defects inside and carries the symmetry of 
' the director and it enables to determine qualitatively the type of the interaction potential. It was first shown there 
that axial symmetry breaking leads to the origin of elastic charge and coulomb-like interaction. But the coat is not 
quantitatively exactly defined. In 11, authors suggested approach by fixing director field on the surface of the 
, imaginary sphere. There multiple expansion of the potential was obtained as well but coefficients differ significantly 
from the results of Q. For instance dipole-dipole interaction was found to be three times weaker, and quadrupole five 
times weaker than in On the other hand authors of [f| recently have measured experimentally both interactions 
and found that results are in accordance with Lubensky at all prediction with about 10% accuracy. This allows to 
confirm Tightness of the main assumptions of [7| for spherical particles though it do not work for particles with broken 
axial symmetry. In this paper we suggest to generalize that results for the case of arbitrary particle's orientation. We 
. , propose new free energy functional that describes interaction of colloidal particles with director field and turns into 
functional of Q for axial particles. It is valid for both strong anchoring case with the presence of topological defects 
and for weak anchoring as well. We show that axial symmetry breaking leads to the violation of SpQ^ — condition 
for spur of the quadrupole moment of particles in contrast to the classical electrostatics. This leads to the principal 
new multipole terms in elastic interaction connected with SpQ^ which are absent in the classical electrostatics. We 
show that the reason of such difference leads in nonlinearities of director field on close distances and in surface driven 
effects associated with anchoring. 

In order to explain subsequent speculations we consider first the case of homogeneous liquid crystal with uniform 
director no and one particle immersed into it. Anchoring of the liquid crystal with the surface of the particle produces 
deformations of the director field around the particle so that director n(R) varies from point to point. In the one- 
constant approximation the total free energy of the system: particle plus liquid crystal has the form: 

F = — J dV [(divn) 2 + (rotn) 2 ] + j> dSW{v • n) 2 (1) 

where W is anchoring strength coefficient, v is the unit normal vector, integration § dS is carried out on the surface 
of the particle, K is the Frank elastic constant. 
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Far from the particle director field variations are small n(R) w (n x , n y , 1), |n M | <l(/j = x, y) and bulk free energy 
has the form: 

F b , Hnear = y /" {(V»( x ) 2 + (Vn a ) 2 } (2) 

which brings Euler-Lagrange equation of Laplace type: 

An, = (3) 
At large distances R in general case it can be expanded in multiples, 

with m q = R a /R; p^u = p^u Q , u : : u = Q""u a Up, SpQ^ — Q^ a and /i = (x,y). This is the most general 
expression for the director field. It differs from usual electrodynamics with —SpQ^/R 3 term. We shall show below 
that it is not zero when the axial symmetry breaking takes place. We note that the multiple expression does not 
depend on the anchoring strength. It is valid on far distances for any anchoring, weak and strong, without topological 
defects or with them. Of course in order to find multiple coefficients we need to solve the problem in the near nonlinear 
area either with computer simulation or Ansatz functions. Let's imagine that we have found all multiple coefficients 
for the particular particle (for instance with help of computer simulation): 20 different values, with 10 for every n^. 
So for every p = (x,y) we have one elastic charge q^, three components of the vector dipole moment p^ and six 
components of the quadrupole tensor Q M (it is symmetric Q°f = Q^ a , but SpQ^ ^ in general ) . How particles will 
elastically interact if we know all these 20 multiple coefficients. We want to obtain analytical interaction potential. 
Our task is to build effective harmonic free energy functional which includes all multiple coefficients and nothing else. 
It must be invariant and it's Euler-Lagrange equations for one particle must have solution in the form of multiple 
expansion (j4j. It is easy to check that such functional has the following form: 

F = K J d 3 x ji^L - 4^(x)rv - 47rj£(x)0 Q rc M - 4 7 rQ^(x)d a d /3 n M } (5) 

In the case of N particles q M (x) = I^g^,i5(x-Xi);^(x) = J2iP°,i S ( x ~ x «); Q"^( x ) = J2i Q^( x-X i) witn i bein S 
the number of the particle. Euler-Lagrange equations for functional ([5]) are Poisson like: 

An, = -47nz M (x) + 4n[d a (pfa)) - d a d p Qf(x)} (6) 

Substitution of this equation into ([5]) leads to the following expression for the total energy of the system: particles 
and liquid crystal: 

F total = -2ttK J d 3 x { 9 „(x)n« + (p„(x) • V)n» + (Q M (x) : V : V)n~} (7) 
where n e ® -exact solution of the ^ . For unlimited space the solution of equation ^ has standard expression: 

r V (x) = J d'x'j^^ k(x') ~ 5>«(x')) + d' a d' Qf( X ')] (8) 

For one particle we have <? M (x) = g M (5(x),p"(x) = p"<J(x), Q""(x) = Q"f 3 S(x) which leads us immediately to the 
general expression 0$. For axially symmetric particles = and the director take the form n x = p-j^ + 6Q^,n y — 
Pr? + QQt§:- In Lubensky notation Q Q = §. 

Actually in the construction of the functional ([5]) we follow Lubensky ideology but with one significant difference. 
In the paper [t| authors tried to build invariants of vector P, tensor C\j and director vector n, which has three 
components, not two. So obtained invariants are of the type — P • n(Vn), (Vn)n • V '(mCijUj) and include nonlinear 
terms of n^. In linear approximation authors find the following main elastic functional: 

F = K [d 3 x!^^- 4vrP z (x)9 A1 n Al - iirC^d^nA (9) 
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This functional is axially symmetric, it correctly describes axially symmetric particles. But in the case of axial 
symmetry breaking it does not valid. In the case of axial symmetry breaking we need to take into account coulomb 
like expression and to introduce two different vectors p^p^ and two tensors Q x ,Q y for each component /i. Using 
them we find new functional |(SJ) which describes system of any particles immersed into NLC. Let's show how ([5]) 
proceeds to the © for axially symmetric particles. In this case p^, = (p,0,0);p = (Q,p, 0) and quadrupole tensors 
are of the form: 

/ Q \ / \ 

(10) 

Then automatically © proceeds to the © with P z (x) = p(x);C zz (x) = 2Q(x). 

It is necessary to note that functional ((SJ) defines energy of the colloidal particle in the director deformations created 
by any external sources. So if we have some deformed NLC the particle will move in the way to minimize ([5]) . 

Director field from N immersed particles will be superposition of deformations J3J from all particles « M (x) = 

E; r^Tf + E* P ^"' ) + E* ^^Ijfc^ - E, f^f • Substituting this director into © and neglecting with 
self energy divergent parts we come to the fact that total energy of the system may be presented as the sum of pair 
interactions: 



- r 





Q\ 







°\ 


Qx= o 





)Qy = 


: 





Q 


\Q 





0/ 




Q 


0/ 



F total = Y,Y, ut1 ^-^) ( n ) 

i j>i 

ipi = U qq + U qd + U dd + U qQ + U dQ + U QQ (12) 



U qq = -4nK^^i (13) 



= _ 47rX ^(PM,r^)-gM,-(p„r^) 

^ = _ 47rj ,(P„-P,,)-3(P r -u,)( PAtJ .u iJ ) (i5) 

u \2ttK I < hb^ii '■ Qm : + : : u »j) (q»,iSpQ»,j + g^jSpQ^ 1 

1 ''<'/ 3r ij J 

cn r ,( u u : 4m : u y)(P M ' u y) ~ ( u 'y : Qm : u «)(P ft i ■ u *i) 

— 60nK — 2 — 

r 

C^QQ = — I'Z'kK < :: Q^'J ~ 20 u ij | Q^.j ■ : u ij + 35(iijj : : Ujj)(ujj : Q^j : Ujj) \ _ (18) 

T ij J 

1 ij 
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where Ujj = |*' _* J | and r,j = |x,- — Xj|. These expressions are the most general elastic pair interaction potentials 
between ordinary particles immersed into NLC. They are valid for particles with topological defects and without 
them, for any type and strength of the anchoring coefficient. Coulomb-like elastic attraction and repulsion between 
inclined ellipsoidal particles was found in [l6j], pjj with help of computer simulation. We consider that Coulomb- 
like elastic interaction takes place between any non-symmetrical particles with broken axial symmetry (for instance 
non-symmetrical triangles or cubes with different anchoring on sides) even in the case of zero torque applied to the 
particles. 

If particles are the same then U q d = 0, UdQ = 0. We see that formulas (fT6 |) -(ri~8 |) contains new terms connected with 
SpQfj. which are absent in the classical electrostatics. Let's analize this result in comparison with electrostatics. In 
electrostatics there is the charge density distribution p(x) in the Poisson equation A(p = — Anp and linear electrostatics 
works inside this density. So inside p(x) there is superposition principle. But in nematostatics there is principle 
difference on small distances - the superposition principle does not work because of nonlinear deformations. So 
in case of strong anchoring we consider that the notion of elastic charge density is unreasonable. That's why in 
nematostatics we introduce the concept of elastic charge q^, dipole moment p M , quadrupole moment Q^f not as 
compound expressions like in electrostatics (full charge is equal the sum of compound charges qtotai — 9*; dipole 
moment as P = '^2 i q%Yi and so on ), but as primary quantities for each particle, the quantities which are just the 
linear asymptotics of exact nonlinear solutions. 

As well in nemato-statics surface effects can make SpQ^ ^ 0. Let's show this on the example of small particles 
with weak anchoring without topological defects. Then deformations are small anywhere n(R) w (n x , n yi 1), |n M | <C 1 
(/i = x,y) and surface energy takes the form §dsW{v ■ n) 2 ps const + 2W <f dsv^v^n^s). Taylor expansion of 
the director concerning center of the particle brings : n^(s) = n^(x a ) + (I ■ V)ro M (a;o) + \{l ■ V) 2 n p (xo) where 
1 is the vector from the center x to the point on the surface s = x + 1. Then we automatically obtain that 
free energy ((T|) takes the form ([5]) with such multiple coefficients: q^ — — -^j? § dsvzv^, p p = — § ds^v^s), 
Qff = — 4^j- § dsvj,v^l a {s)lp{s). These quantities are of course approximate as Taylor extension is not exact but 
we consider that qualitative estimates may be correctly performed. Let's note that substituting them into Uij we 
naturally come to the results of papers 0, [l(| . 

So we may clearly see that SpQ^ = — § dsv^v ^l 2 (s) sa (l 2 ) 7^ 0(s) if qy, 7^ 0, so that SpQ^ and are 
approximately proportional to each other ( here (l 2 ) is average square of the distance from the surface to the center of 
the particle) . As the elastic charge occurs when axial symmetry is broken [l(| we come to the conclusion that in that 
case as well SpQ^ 7^ 0. In the paper authors have found Coulomb-like elastic interaction between two inclined 
ellipsoidal particles using computer simulation. As well there analytical expression for the elastic charge for the case of 
small anchoring was obtained using q^ = — § dsu^v^. Then it was found that q x = — F(b, s) cos sin # cos (9, 
q y = ~-^[^F{b,e) sin tp sin 9 cos 9, where 9 is the angle between long axis and n , ip is the angle between the horizontal 
projection of the particle and direction of axis x, b is the length of the long axis and e is the eccentricity e = yl — c 2 /b 2 
with c being the short axis of the ellipsoid and: 

F(b, e) = Trb 2 Tl - T 2 1 (2s 2 - 3) + (3 - 4e 4 ) I . (19) 

In the case of small eccentricities it has the asymptotic expansion F(b, e) £ ^>o — — -^7r6 2 e 2 + o(e 2 ). Here we come to 
an interesting conclusion, that the elastic charge goes to zero as square of the eccentricity q^ cx e 2 . Then we may 
estimate SpQ^ w q^ (l 2 ), where (l 2 ) is something average between b 2 and c 2 . 

In this paper the new free energy functional which describes elastic interaction between arbitrary colloidal particles 
and NLC has been proposed. It contains multiple coefficients of each colloidal particle and proceeds to Lubensky 
functional for axially symmetric particles. This approach enables to find clastic interaction energy between emersed 
particles in the most general case for arbitrary anchoring strength on the surface. The procedure to find interaction 
energy consists of two steps - first, multiple coefficients for director field are found using different methods like 
computer simulation, variational ansatz and so on, second - these coefficients are used in the obtained analytical 
expressions for elastic potentials. It was first time shown that spur of the quadrupole moment tensor is different from 
zero SpQ^ 7^ in the case of axial symmetry breaking of the director field around colloidal particle in contrast to 
classical electrostatics. 
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